Fractal Weyl laws for amplified states in PT-symmetric resonators 
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We find that in nonhermitian PT-symmetric systems (as realized in resonators with balanced 
absorption and amplification), a mechanism based on quantum-to-classical correspondence reduces 
the occurrence of strongly amplified states. The reduction arises from semiclassically emerging 
hierarchical phase-space structures that are associated with the coupling of the amplifying and 
absorbing regions (forward and backward-trapped sets and their complements), and amounts to a 
generalization of the fractal Weyl law, earlier proposed for ballistically open systems. In the context 
of the recently introduced class of PT-symmetric laser-absorbers, this phenomenon reduces the 
number of states participating in the mode competition. 

PACS numbers: 42.55.-f, 03.65.-w, 05.45.Mt, 42.25.-p 



Nonhermitian systems which can still possess a real 
spectrum have been realized very recently in optical [l], Q 
electronic Q and microwave Q settings. These experi- 
ments utilize the concept of VT symmetry 0, — when 
an amplifying system is coupled symmetrically to an ab- 
sorbing system, this can result in a precise balance of 
loss and gain for individual modes, which then have 
a real energy (the time reversal operation T converts 
loss into gain and vice versa, while the parity opera- 
tion V interchanges the amplifying and absorbing parts 
of the system). For small absorption and amplification 
rates /i, this balance is robust, but if these rates are too 
large, pairs of real eigenenergies bifurcate into complex- 
conjugate pairs. This phenomenon is known as sponta- 
neous PT-symmetry breaking and leads to a range of 
remarkable switching effects [1, i, @-@| . In particular, 
PT-symmetric optical resonators have been predicted to 
form a new class of lasers EMI, 

which show an addi- 
tional spectral peculiarity: lasing modes are degenerate 
with perfectly absorbing modes 11-13|. In the setting of 
these lasers, the focus shifts to the most strongly ampli- 
fied modes, as these are best placed to overcome external 
losses, win the mode competition, and thus determine 
the laser threshold and radiation characteristics. 

Here we report that the most strongly amplified modes 
in these PT-symmetric resonators show the hallmarks of 
yet another distinctive spectral feature — a fractal Weyl 
law, by which the number of these modes is systemati- 
cally reduced when the resonator dimensions are scaled 
up to become much larger than the wave length (the 
semiclassical limit, corresponding to an effective Planck's 
constant h <C 1). Fractal Weyl laws have been origi- 
nally established for passive leaky systems with chaotic 
underlying classical dynamics, where long-living states 
are supported by a fractal repeller — the number of these 
states then scales as h dn , where du < d is the fractal di- 
mension of the repeller in the d-dimensional phase space 
[14l4l8|. In order to uncover the analogous effects in 
the PT-symmetric setting, we establish a quantum-to- 
classical correspondence lil-[2o| of various components 



of the spectrum to specific regions in the classical phase 
space, with phase-space volumes reflecting the propor- 
tions of these components in the spectrum. We find that 
the strongly amplified states are supported by a hierar- 
chical structure associated with the coupling of the am- 
plifying and absorbing regions (the backward-trapped set 
in the amplifying parts of the system, which forms part 
of the classical repeller). The increasing phase-space res- 
olution in the semiclassical limit reveals the sparse fractal 
nature of the backward-trapped set and thus results in 
the reduction of the number of strongly amplified modes 
mentioned above. 

Model for nonhermitian VT- symmetric resonator 
dynamics. — A model of the dynamics of waves in a res- 
onator which displays the required VT symmetry can be 
set up as follows (see Fig. [I]). We divide the Hilbert space 
into two subspaces L, R, each of dimension M, with the 
first subspace representing the absorbing region in a 'left' 
resonator and the second subspace representing the am- 
plifying region in a 'right' resonator. In each resonator, 
the ballistic chaotic dynamics over a finite time interval 
At = 1 is described by an M x M-dimensional unitary 
time-evolution operator F L = F and F R = [F -1 ]* = F T , 
respectively, where the stated relation between Fl and 
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FIG. 1: (a) Sketch of a PT-symmetric resonator, consisting 
of an absorbing resonator (left) coupled to an amplifying res- 
onator (right). The unitary operator F describes the internal 
dynamics in each resonator, the factors e ±fl model amplifica- 
tion and absorption, and the projectors P (onto the interface) 
and Q (onto the walls) describe the coupling, altogether re- 
sulting in the quantum map ©. (b) The Hilbert space (and 
the classical phase space) divides into an absorbing and an 
amplifying part, coupled in a subspace (shaded gray). 
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Fr ensures VT symmetry [To|, EI]. In each time step, 
the wave amplitude in the absorbing region is reduced 
by a factor exp(— /i), while in the amplifying region it 
is enhanced by exp(/i). This breaks the unitarity of the 
time evolution, but because of the matching rates again 
respects VT symmetry. Finally, the ballistic coupling 
between the subspaces is described by a coupling matrix 
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diag(l i: ^,0 1 __0). (1) 

N M-N 



where P projects onto the interface between the two 
resonators (with N channels connecting the resonators), 
while Q = 1m — P represents the projector onto the walls 
of each resonator. 

With these specifications, the time evolution of the 
composed system can be written as 











(2) 



tary matrix \[C — 



which 



where we symmetrized the coupling by means of the uni- 
2~ 1 ' 2 P + Q -i2~ 1 l 2 P 
-i2~^ 2 P 2~ 1 / 2 P + Q 
squares to C. The operator T acts on 2M-dimensional 
vectors ip = Q^), with the M entries in iJjl giving the 
wave amplitude in the absorbing subsystem, and the re- 
maining M entries in ipR giving the wave amplitude in 
the amplifying subsystem. 

The spectrum of the system is now obtained from the 
eigenvalue equation 



(3) 



where E n = e n — iT n /2 are quasienergies. The VT sym- 
metry of the quantum map is manifest by the relation 



(4) 



where the parity operator is of the explicit form V = 

1m \ _ ^ M ^ a ^ p interchanges the am- 

1m U J 

plitudes of the subsystems). The spectral properties as- 
sociated with VT symmetry are now embodied in the 
characteristic polynomial s (A) = det ( T — A), of degree 
2M. Due to VT symmetry, this polynomial exhibits the 
mathematical property of self-inversiveness (2ll |22|: 



s(l/\*) = [\-' M s(\)]*s(0), 



(5) 



where 5(0) = det J 7 = (detF) 2 is a unimodular complex 
number. The eigenvalues are the roots of the secular 
equation s (A) = 0. For each eigenvalue A n , we are thus 
guaranteed to find the eigenvalue An = [A" 1 ]*. It is pos- 
sible that n = n; then |A n | = 1, which means that the 
quasienergy E n = iln\ n is real. For n ^ n, we obtain 
complex quasienergies, which thus still appear in pairs 



(E n = E*). If fx = 0, T is unitary, and all quasiener- 
gies are real. In the semiclassical limit of large M = h~ x 
and fixed classical coupling strength N/M, and with \i 
increasing from zero, one expects a quasi- monotonous in- 
crease of the fraction of complex energies. 

Dynamical signatures of spontaneous VT-symmetry 
breaking. — We now investigate how the transition from 
a real to a complex spectrum depends on the dynamics of 
the system. In the resonator setting, nonintegrable dy- 
namics arise from resonator walls of a generic shape. In 
the quantum map ([2|), this can be modeled by choosing 
a suitable operator F for the internal dynamics. In an 
optical setting, the passive system is time-reversal invari- 
ant, which furthermore dictates F — F T . We implement 
these features via a paradigm of wave dynamics with a 
nonintegrable classical limit, the kicked rotator [23I, l24j. 
The time-evolution operator of this system is given by 



-+COS - 



(6) 

where the kicking strength k controls the dynamics from 
classically integrable (k = 0) to globally chaotic (k > 7); 
the classical map is q r = q+p+ (k/An) sin(27rg) (mod 1), 
p' = p + (/c/47r)[sin(27rg) + sin(27rg / )] (mod 1). We focus 
on the chaotic cases with k — 8, fix the coupling of the 
absorbing and amplifying regions by setting the inverse 
average dwell time in the amplifying or absorbing regions 
(the Thouless energy) to Et = N/M = 1/5, and analyze 
the spectrum (obtained by diagonalizing T) as a function 
of the amplification rate \i and the system size M. 

We start with the distribution of decay rates, encoded 
in Imi^, where large positive values indicate strong am- 
plification and large negative values indicate strong de- 
cay. The left panels in Fig. [2] show histograms of this 
quantity for three representative values of \i. For compar- 
ison, the right panels show results from random-matrix 
theory, with F taken from the circular orthogonal en- 
semble [25[- All histograms are symmetric because com- 
plex energies appear in conjugate pairs, as imposed by 
VT symmetry. Furthermore, the histograms display a 
sharp peak at IvaE = 0, which arises from the states 
with real energies. This peak decreases with increas- 
ing systems size M, indicating that the proportion of 
such states vanishes in the semiclassical limit. As re- 
gards to this peak, there is good correspondence between 
the dynamical system and random-matrix theory, which 
predicts that the spectrum essentially turns complex at 
fi c = VN/M <^E T [It 

For jjL = O.bEr (top panel), this agreement between 
both models also extends to finite values of ImE, even 
though some difference are noticeable. The differences 
become more marked as \i exceeds Et (middle and bot- 
tom panels) — in particular, in comparison to random- 
matrix theory, the dynamical system possesses a signifi- 
cantly reduced number of strongly amplified states, with 
large YvaE. The insets in the bottom panel focus on the 
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FIG. 2: (Color online) Distributions of ImE in a "PT- 
symmetric resonator with Thouless energy Et = N/M = 1/5 
for various amplification rates \i and system sizes M, nu- 
merically obtained from Eqs. J2]), J3} with internal dynamics 
modeled via the kicked rotator ((6} with k = 8 (left) and via 
random-matrix theory (right). Significant differences occur 
when \i exceeds Et- For the kicked rotator, the proportion 
of strongly amplified states decreases systematically with in- 
creasing M, as highlighted in the inset of the bottom panel. 



peak associated with these states. In the dynamical sys- 
tem, their proportion decreases systematically as M in- 
creases, while in random-matrix theory their proportion 
does not change. 

In order to quantify this systematic reduction we con- 
sider the fraction of states with ImE > /i/2, denoted as 
/> and plotted as a function of \i in the left panel of Fig. 
[3l In random-matrix theory, this curve becomes indepen- 
dent of M when M is large. For the dynamical system, 
however, />(/i) drops as M increases, even at the largest 
computationally accessible system sizes [M = 8 000, cor- 
responding to a matrix dimension 16 000 for J 7 ). As 
shown in the right panel this drop follows a power-law 
/>(/x) oc M~ a , with a = 0.102(9) (/i = E T ), a = 0.079(6) 
\li = 2E T ), and a = 0.069(8) (/i = 2E T ). This data con- 
firms that deviations from random-matrix theory set in 
as /i exceeds Et- That the deviations increase with M 
points to a semiclassical origin, which we identify next. 

Following previous investigations of general non- 
hermitian systems that addressed the issue of state 
nonorthogonality 



18, 



[I9j, we identify the semiclassical 
support of the eigenstates in different components of the 
spectrum by constructing an orthonormal basis (j) n of 
the subspace spanned by the states, and then summing 




FIG. 3: (Color online) Left panel: Fraction />(/x) of strongly 
amplified states with Imi^ > fi/2, as a function of \i. The 
random-matrix curve (RMT) is evaluated at M — 8000, while 
the curves for the dynamical system (specified in Fig. [2} cor- 
respond to M = 400, 1000, 2000, 4000 and 8000. The double- 
logarithmic plot in the right panel demonstrates power-law 
scaling /> oc M~ a at fixed \i — Et (cl = 0.102(9)), 2 E T 
(a = 0.079(6)), and 5 E T (a = 0.069(8)). 

(a) absorbing region amplifying region (b) backward-coupled 




FIG. 4: (Color online) (a) Husimi phase-space distributions 
of the support of strongly amplified states (ImE > /x/2, 
'neutral' states (llmi^l < fi/2, H n ), and strongly decaying 
states (Imi? < — /x/2, obtained for the kicked-rotator 

model of a T^T-symmetric resonator with k = 8, M = 1000, 
Et = 1/5, [i — 2Et. (b) Classical phase-space regions which 
are coupled to the interface by 1 — 4 steps, backwards in time 
(top panel) or forwards in time (bottom panel). 



the Husimi phase-space representations H(q,p) of the 
basis states, defined by the squared overlap |(^,p|^> n )| 2 
with minimal-uncertainty wavepackets centered at phase 
space position (q,p). Figure a) shows these distribu- 
tions for /jl = 2 Et and M = 1000. The top panel depicts 
the support H > of the states with ImE > /i/2, which 
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contribute to the quantity />(/i) discussed above. The 
middle panel shows the support H n of 'neutral' states 
with |Imi£| < /i/2, while the bottom panel shows the 
support H < of quickly decaying states with ImE < 
—/i/2. The distributions H > and H < are related by VT 
symmetry, which interchanges the absorbing and ampli- 
fying regions and maps the momentum p — >• — p; H n is 
invariant under these operations. 

The systematic structures observed in these distribu- 
tions have a specific classical dynamical interpretation, 
which follows from the plots shown in panel (b) of Fig. 
IU The top panel in (b) shows backward- coupled regions 
which are populated from the interface after 1—4 times 
steps of the classical dynamics. Including higher itera- 
tions, this hierarchical structure defines the complement 
of the backward-trapped set, a fractal which contains all 
points whose backwards-in-time dynamics does not meet 
the interface. The lower panel in (b) shows the analogous 
forward- coupled regions, which propagate towards the in- 
terface and define the complement of the forward-trapped 
set (the time reverse of the backward-trapped set). 

Coming back to panel (a), we now recognize that in 
the amplifying part of phase space, the strongly ampli- 
fied states mainly populate the backward-trapped set, 
and thus concentrated in regions of maximal dwell time. 
In the absorbing part of the system, these states are sup- 
ported by the backward-coupled regions, which signifies 
minimal dwell time in this part of the system. A power- 
law suppression ex M~ a of the number of these states now 
follows from the fractal nature of the backward-trapped 
set, which is resolved in more detail in the semiclassical 
limit (according to the shrinking size 1/M of a Planck 
cell in phase space). The relevant fractal repeller dimen- 
sion can then be estimated from the data in Fig. [3] as 
dn = 2 — a (the precise definition of this dimension in 
the context of fractal Weyl laws is an open problem). 
For larger values of /i, /> picks up larger parts of the 
bulk of the distribution P(lmE) : thus reducing the slope 
and overestimating da- We note that there are quan- 
tum fluctuations as a function of M, and classical finite- 
size effects as the dwell time E^ 1 = 5 is not asymptot- 
ically large. For the strongly decaying states, the same 
argumentation carries over to the forward-trapped set. 
Furthermore, the neutral states are supported by the 
backward-coupled regions in the amplifying part, and the 
forward-coupled regions in the absorbing parts, thus re- 
sulting in an (approximate) balance of amplification and 
absorption which becomes exact for states with Im E = 0. 

In conclusion, the spectral properties of T^T-symmetric 
systems with complex wave dynamics show clear signa- 
tures of the underlying classical dynamics, and in par- 
ticular, the details of the coupling between amplifying 
and absorbing regions. Strongly amplified states are sup- 
ported by the backward-trapped set in the amplifying 
region, while states with real or almost real energies are 
supported by regions that are well connected to both the 



amplifying and the absorbing parts of the system. The 
fractal nature of the backward-trapped set results in a 
systematic reduction of strongly amplified states, which 
becomes more marked for large system sizes, and follows 
the characteristic power-law dependence of a fractal Weyl 
law (earlier predicted only for leaky systems with ballistic 
openings). 

The experimental observation of the ordinary fractal 
Weyl law has proven a challenge as it addresses decay- 
ing quasi-bound states. In contrast, the fractal Weyl 
law uncovered here applies to amplified and possibly las- 
ing states. This phenomenon should generally affect the 
properties of the discussed class of T^T-symmetric laser- 
absorbers, where the strongly amplified states dominate 
the mode competition. However, conceptually our find- 
ings also carry over to other classes of nonhermitian dy- 
namical systems, including passive realizations of VT 
symmetry [H, H| , as we U as systems which combine am- 
plification and absorption in a non-symmetrical fashion 
[27j. The main requirements are multiple scattering on 
long time scales and ballistic wave propagation on short 
time scales, including a ballistic coupling between the re- 
gions of different amplification or absorption rate; the 
most amplified (or least decaying) states are then semi- 
classically supported by the backward-trapped set in the 
most amplifying (or least absorbing) regions. 
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